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Abstract. 

It is shown that, under suitable conditions, involving in particular the existence of analytic 
constants of motion, the presence of Lie point symmetries can ensure the convergence of 
the transformation taking a vector field (or dynamical system) into normal form. 
A short version of this paper has been submitted, as a Letter, to the Jour-§ 
nal of Phys . A. 



1. Introduction. 

The technique of transforming a system of first-order ordinary differential equations 
(also called dynamical system) into "normal form" is an old and well known method 
of investigation, going back to the classical work of Poincare (and then of Dulac and 
Birkhoff), and developed in more recent times by several authors (see e.g. [2,3,5-7,22] and 
ref. therein). Its connection with symmetry properties (precisely. Lie-point symmetries 
(see [20,21])) of the dynamical system have been also pointed out [1,6,11,12,16]. One of 
the main troubles with this procedure is given by the problem of the convergence of the 
normalizing transformation: it is known in fact that these transformations are performed by 
means of recursive techniques, and only special conditions can ensure their convergence [7] . 

A possible way-out which is usually adopted is that of considering these transforma- 
tions "up to some finite order", i.e. of considering "approximate normal forms" (and, 
correspondingly, approximate solutions, approximate symmetries, etc. of the given sys- 
tem; see e.g. [4,14]). In this paper, instead, we are dealing with the case of converging 
normalizing transformations, and with some conditions ensuring convergence. Our discus- 
sion is related to a theorem by Markhashov [19]: even if the original proof of this theorem 
has appeared to be not complete, a similar result has been proved at least in some partic- 
ular case (see [8]). A very remarkable result in this context is given in a recent paper by 
Bruno and Walcher [10], in the case of 2-dimensional systems. The main idea is that the 
presence of some symmetry of the dynamical system can ensure - under suitable conditions 
- the convergence of the normalizing transformations. 

Our paper presents some considerations in the same direction. In particular, we give 
a direct proof of a "Markhashov- type" theorem in a well definite and simple case, and 
analyze some applications of this idea, which include a generalization of the Bruno and 
Walcher result [10]. 



2 



2. A " Mar khashov- type" theorem. 



We will consider dynamical systems (DS) of the form 

ii = f{u) u = u{t) e i?" (1) 

where ii = dw/dt, / is assumed to be analytic in a neighbourhood of it = 0, with /(O) = 0, 
and 

f{u) = Au + F{u) (2) 

where the matrix A = (V/)(0) is assumed to be nonzero and diagonalizable. Let us remark 
that most of the results below could be extended to the non-diagonalizable case, apart from 
some complications in the notations and statements (see [1,13]). 

As well known [2,7], a normalizing transformation is a nonlinear formal transformation: 

u ^ u = u + . . . (3) 

transforming (1-2) into a new DS which we write in the form (to avoid cumbersome nota- 
tions, we will denote by u both the "original" and the transformed coordinates) 

u = f{u) ^Au + F{u) (4) 

where the nonlinear part F{u) is in "normal form" (NF). To define this notion, one in- 
troduces in the space of analytic functions, defined in a neighbourhood of tt = 0, the 
Lie-Poisson bracket 

{/, 9}k = (/ ■ V)^fe -{g-V)fk (/c = 1, . . . , n) (5) 

and, given any n x n matrix A, the "homological operator" A 

A{f) = {Au, /} = {Au) -Vf-Af. (6) 

Then, a nonlinear vector function F is said to be in NF with respect to A (or resonant 
with A) if 

AiF) = (7) 

In the basis where A is diagonal, with eigenvalues ai,...,a„, a monomial Fk{u) = 
u^^ ■ . . . ■ u^" of degree j (with integer numbers such that rrii = j, rrii > 0) is 
resonant if (m, a) = ^j^rriiai = ak- As well known, the relevance of the above definitions 
is essentially due to the fact that, given a vector function /, all nonresonant terms can be 
removed by means of a formal coordinate transformation. 
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^^/■Vp-O. (10) 



We also say that a vector function 

g{u) ^ Bu + G{u) (8) 

is a (Lie-point time-independent) symmetry for the DS (1-2) if 

{/,^} = 0. (9) 

In terms of Lie algebras, one says that the vector field operator g ■ V generates a symmetry 
of the DS. 

A scalar function p = p{u) is a (time- independent) constant of motion (or first integral) 
for the DS (1-2) if the Lie derivative along / vanishes: 

dp 
dt 

The above definitions of symmetry and of constant of motion can be clearly applied both 
to analytic functions and to formal power series. 

Our discussion needs few preliminary results, some of which are rather simple or well 
known; however, for clarity and completeness, we give all of them: some of these introduc- 
tory results may also have an independent interest. 

Lemma 1. If is a symmetry for the DS, and p a constant of motion for it, then also 

h = p g 

is a symmetry for the DS. More precisely, the algebra of the symmetries of a DS is module 
over the constants of motion of the DS. 

The first part of this Lemma is immediate; the other statement describes the general 
property [13,15,16,21] of the solutions g to the system of PDE's (9) which gives the sym- 
metries of (1-2). See later on for some remark concerning the number of the "admissible" 
constants of motion. A fundamental step in the discussion is provided by the following 
Lemma. 

Lemma 2 [1,12,23]. If the DS (1-2) admits a symmetry g (8), and / is put in NF 

f ^ f = Au + F, FGKer^ (11) 

by a formal normalizing transformation, then g is transformed into a new form (not nec- 
essarily normal and possibly formal) which we denote by (we reserve the notation ^ only 
to NF) 
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g^g = Bu + G (12) 

satisfying (together with {/, g} — 0, foUowing directly from (9)) 

G e KevA or equivalently {Au, G} = Q (13) 

F e Ker^B or equivalently {Bu, F} = (14) 

where B is the homological operator B{-) — {Bu, ■}. 
Proof (a sketch). Expanding F in formal power series: 



J>2 



where Fj are homogeneous polynomials of degree j, and expanding in a similar way G, 
one has from {/, ^r} = at order 1 

[A,B]^0 (15) 

and at order 2 

{Au,G2} = {Bu,F2} (16) 

Applying to this equality the operator A, one has A'^{G2) = which gives, thanks to the 
assumption on the matrix A (see also [12]), A{G2) — 0, and from (16), also B{F2) = 0. 
Iterating the procedure, one obtains (13-14). 

Remark [12]. By means of a further formal transformation, one may also normalize g 
thus obtaining the "joint" NF: 

F, G e Ker^ n KerB . (17) 



An immediate but very important consequence of the notion of NF and of Lemma 2 
(in particular eq. (14)) is the following proposition. 

Proposition 1. If the DS is in NF (11), then it admits the linear symmetry Au. If the 
original DS (1-2) admits a symmetry g = Bu + G, then the normalized DS (11) also admits 
the linear symmetry Bu . 



Let us now recall the two basic conditions which ensure that a vector function / = 
Au + F can be put in NF by a converging transformation [7] , namely: 
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Condition "A": there is a coordinate transformation changing / to /, where / has the 
form 

/ = Au + a{u)Au 
and a{u) is some scalar- valued power series (with a{0) = 0). 

Condition "oj": let ojk = min Kg-, a)\ for all positive integers Qi such that Yl7=i ^* ^ 
{q, a) 7^ 0: then the series 

oo 

^2-'=lnc^fc 
k=l 

is convergent. 

The first condition will play a key role in our discussion; the other one is a weaker condition, 
controlling the appearance of small divisors [7] , and we explicitly assume that it is always 
satisfied here. In particular, it is satisfied in the cases considered in sect. 4. 

We are now in position to give a simple and direct proof of the following result. Let 
us preliminarily note that any DS admits an obvious symmetry, namely g = f, which is in 
fact the generator of the dynamical flow. Accordingly, when considering the symmetries 
of a DS, it is always understood that none of them (nor their linear combinations) is 
proportional to /. Also, we may clearly exclude the case (which may be considered here as 
"trivial") that the DS, once in NF, takes the form u = (l + a{u))Au: in this case indeed 
the Condition "A" is satisfied, and the convergence of the normalizing transformation 
automatically guaranteed. 

Theorem 1. Assume that the DS (1-2) possesses a finite number £ (> 1) of analytic 
symmetries gj = BjU + Gj, where all the matrices Bj are linearly independent (and not 
zero). Assume also that, once in NF, the DS admits exactly £ linearly independent (possibly 
formal) symmetries. Then, there exists a converging normalizing transformation for /. 
Proof. Let us start by writing the DS into NF by means of some formal normalizing 
transformation 

f^f = Au + F 
Under this transformation the £ symmetries become 

9j 9j = BjU + Gj (18) 
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According to Prop. 1, the linear parts BjU together with Au are {£ + 1) symmetries of 
the DS in NF; excluding the trivial case / = Au, this implies that there must be a linear 
combination of the Bj which is equal to A: it is not restrictive to assume e.g. Bi = A; 
consider then 

gi=Au + Gi 

By assumption, our problem in NF admits exactly £ symmetries (plus the trivial one /), 
therefore gi must be a linear combination (with constant coefficients) of the symmetries 
at disposal, i.e. of sj = Bju and /. But the linear part of gi is just Au, and this forces 
Gi = 0, or 

gi = Au 

This implies that Condition "A" is satisfied by the transformation of the symmetry gi = 
Au + Gi — > 5^1 = Au; then there is a normalizing transformation which is convergent. 
Under this transformation, / is transformed into / = Au + F which is in NF, according to 
Lemma 2 (indeed F e Ker^, being gi = Au + Gi). 

The above theorem looks quite "formal" and not easily applicable in concrete cases: in 
fact, it may be difficult to check in practice that the required properties of the NF (which is 
usually "a priori" not explicitly known) are verified. In the two next sections, we will give 
more concrete versions of this result and study some cases in which the above hypotheses 
can be fulfilled. 

3. The number of constants of motion and symmetries of a DS. 

First of all, let us remark that one of the crucial hypotheses of Theorem 1 is that the 
DS in NF admits a finite number of independent symmetries. According to Lemma 1, 
the finiteness of this number depends in an essential way on the number of independent 
constants of motion of the DS. In fact, it is clear that the presence of constants of motion of 
the problem in NF precludes the application of our above argument: indeed, the nonlinear 
part Gi of the symmetry gi could be obtained as a combination of / and of the other Bju 
multiplied by suitable constants of motion, and Condition "A" for the symmetry gi fails 
to be verified. 

To carefully discuss this point, let us recall the two following relevant results. 
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Lemma 3 [10,23]. If / = Au + F is in NF, i.e. F e Ker^, then any (formal) constant of 
motion p of the DS = / is also a (formal) constant of motion of the linear part ii = Au. 

Lemma 4 [13,16]. Given the matrix A, the set Ker^ of terms F{u) resonant with A, is 
given by M{iJ,{u))u, where M is the most general matrix such that [M, ^] = and its 
entries My are functions of the constants of motion /i = fi{u) of the linear system u = Au; 
or also, choosing a basis Mj in the space of the matrices M commuting with A, the most 
general F in NF is 



with fij{0) = 0. 

According to Lemma 3, there are two essentially different ways in which a DS may 
admit a finite number of constants of motion, namely: 

1) the linear part u = Au admits a finite number of linearly independent constants of 
motion (then the same is true, a fortiori, for the full DS); 

2) the linear part does admit infinite constants of motion (functionally dependent, of 
course), but only a finite number of them (possibly none) is admitted by the full DS. 

Both these cases deserve some remarks. 

Consider, as an example of case 1), a 2-dimensional system with A — diag(l, 2): then 
a constant of motion is p = u\/u2. This is not an analytic nor formal constant of motion, 
but it can be admitted in this context, because we are interested in analytic symmetries, 
which are described by vector functions of the form g = Bu + G, i.e. with ^(0) = 0: then 
pg may be analytic even if p is not. In the hypotheses of Theorem 1 we assumed Bj 7^ 0, 
then the only admitted constants of motion are of the type u^^ • . . . ■ u^'^ where one at 
most of the may be —1. Therefore, only a finite number of independent constants of 
motion are admitted in this case. A situation of this type happens e.g. when the resonant 
eigenvalues are real and have the same sign. From the point of view of the convergence 
of the normalizing transformations, this case is actually more conveniently treated with 
the notion of the Poincare domain (which provides direct criteria of convergence [2]); the 
present discussion offers simply a different approach to the problem. 

A more interesting situation occurs in case 2), which we are going to consider. 

Assume e.g. that the DS takes in NF the following special form (where — in the sum 
appearing in (19) — only the matrix A and some other matrix M, commuting with A and 




(19) 



not proportional to A, are present): 



u = f = Au + a{u)Au + ijl{u)Mu 



(20) 



According to Lemma 3, any constant of motion of (20) is also constant of motion of the 
linear problem ii = Au: therefore, to find a constant of motion of (20), we start with a 
constant of motion p of the linear problem, i.e. Au ■ Vp = 0, then one gets 



Now, Mu ■ Vp = is verified if and only if p is a constant of motion of the linear problem 
ii = Mu; therefore, if there are no common constants of motion of the problems ii = Au 
and ii = Mu, one has dp/dt|^7^ 0. In particular, this happens if M = / (= Identity) and 
the eigenvalues of A are nondegenerate: indeed the constants of motion of ii = lu arc of 
the form Ui/uj (or functions thereof), but the assumption on the eigenvalues of A excludes 
that these can be constants of motion of ?i = Au. Then we can state the following: 

Theorem 1'. Assume that: i) the DS (1-2) possesses a finite number £ (> 1) of analytic 
symmetries gj = Bju + Gj, where all the matrices Bj are linearly independent (and not 
zero), and where £ is precisely the number of the linearly independent linear symmetries 
admitted by the DS once in NF; 

ii) once in NF, the DS has the form (20), and the two linear problems ii = Au and ii = Mu 
do not admit common constants of motion. 

Then, the DS can be put in NF by means of a convergent normalizing transformation. 

Theorem 1". Condition i) in Theorem 1' can be substituted by one of the two following: 

11) the DS (1-2) admits a nontrivial analytic symmetry g = Bu + G such that 5 7^ is 
proportional to A; 

12) the DS (1-2) admits an analytic symmetry g = G{u) with vanishing linear part and G 
not proportional to F. 

Remarks, a) Assumption ii) in Theorem V ensures that the NF (20) has no other 
symmetry apart from the linear ones BjU. Instead, in the "trivial" case that Condition 
"A" is satisfied by the DS, so that the NF has the form ii = a{u))Au, any constant of 
motion of the linear problem ii = Au is clearly also constant of motion of the full problem 
in NF. Then, in this case — in contrast with the case covered by Theorem 1' — the DS 



dp 
dt 7 



f ■Vp = ii{u)Mu ■ Vp 



(21) 
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in NF does admit, in general, an infinite number of (functionally dependent) polynomial 
constants of motion, and of symmetries as well. 

b) In the assumption ii) of Theorem 1' the request that no common constants of motion 
are present can be substituted by the request that the only common constants of motion 
are rational functions of degree 0. Our final example is an application of this assumption. 

c) In the case the DS satisfies condition {2) of Theorem 1", it is clearly sufficient to consider 
the new nontrivial symmetry g = G + f = Au + F + G to recover ii). The argument then 
proceeds as in Theorem 1. 

4. Applications and examples. 

We want first to show that the theorem by Bruno and Walcher [10] for 2-dimensional 
DS can be reobtained as a corollary of the above approach. We have in fact: 

Corollary. Consider a 2-dimensional DS ii = Au + F such that the eigenvalues ai, 02 of 
the matrix A satisfy a relation kiai + A;2«2 = where fci, /c2 are non-negative relatively 
prime integers, not simultaneously zero. Then, if the DS possesses an analytic symmetry, 
it can be put in NF by means of a converging transformation. 

Proof. Let ii = Au + F the DS, and g — Bu + G its symmetry. If i? = then consider the 
new symmetry g + f = Au + (F + G), so that the linear part is now ^ 0. From [A, M] = 
and observing that the eigenvalues are in this case necessarily distinct, one has that M 
must be a combination of A and /. Then, thanks to Lemma 4, the DS in NF takes just 
the form of (20): 

u = Au + a{u)Au + ii{u)u 

If /U = 0, the NF satisfies Condition "A" and there is a convergent normalizing transfor- 
mation. If |U 7^ 0, the same result follows from Theorem 1': indeed, there are no constants 
of motion for the NF, and the only admitted symmetry is Au; then the argument proceeds 
just as before. 

It can be remarked that in the case of dimension n = 2, one gets directly i = 1, and 
the occurrence of the special form (20) of the NF is guaranteed by Lemma 4. Then, all 
assumptions are automatically fulfilled in dimension n = 2. Instead, if n > 2 and the NF 
has the general form (19), it is clear that a constant of motion p of the linear part ii = Au 
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may be a constant of motion of the full problem u = Au + F, even if it is not a constant 
of motion of each one of the single problems u = MjU. 

One of the possibilities which can guarantee the special form (20) of the NF, and allow 
us to repeat the above argument on the number of constants of motion, is the presence 
of some additional symmetries gj = BjU + Gj: thanks to Lemma 2, one has that F must 
satisfy {BjU, F} = 0, and this condition may exclude some of the matrices Mj in the 
expression (19). The next example will show this possibility, and will also be a good 
illustration of the above discussion. 
Example. 

Consider the space i?^"^, and put u = {xi, . . Xm, J/i, • • • , J/m) £ -R^"^; assume that a 
Lie group F acts "diagonally" on both the m-dimensional spaces of the vectors x and y 
through the same linear representation V: i.e. x ^ x' = Vx, y ^ y' = T>y. Consider then 
a DS of the following form 

u = Au + F (22) 

where 

and assume that F{u) admits the symmetries -BjW, where Bi are the (matrix representatives 
in the direct sum P © "D of the) Lie generators of this group F (the linear part Au fulfils 
this symmetry requirement, so that the full DS (22) admits this symmetry). Let us assume 
here for concreteness (the general case could be relevant for the study of Hamiltonian DS, 
see [17,18], but we do not consider here this situation) the case m = 3, F = SO(2t) and V 
its fundamental representation, and consider the DS 

u = Au-\-p{u)u (23) 

where p{u) is an analytic function depending on the quantities, thanks to the 5'0(3) 
symmetry, x"^ — {x,x), y^ = {y,y), x ■ y = {x,y) (the parentheses stand for the scalar 
product in R^). Once in NF, this DS takes the form 

u = Au + aAu + iiu (24) 

where a and are functions of = + only, thanks to Proposition 1, which ensures 
that the linear symmetry Biu is preserved, and to Lemma 4 as well. This NF has precisely 
the special form of (20), and it is easy to check that there are no constants of motion 
for (24), apart from 0-degree rational functions as {x\ + yi)/{x2 + J/i)- 1^ is important 
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to notice here that, if our problem would not possess the symmetry SO (3), the NF (24) 
would contain many other matrices Mj (according to Lemma 4), but that it is precisely the 
presence of the symmetry SO{S) which forces the NF to contain only A and the identity. 
This seems to confirm the conjecture [9,10,17,18] that the presence of a "sufficient" number 
of symmetries may be an essential request in order to guarantee the convergence of a 
normalizing transformation. Let us now assume (this example can in fact be viewed as a 
multi-dimensional extension of an example given in [10]) that in the original DS (22) the 
function p is a homogeneous polynomial of degree k of the quantities x'^, y^, x ■ y: then 
the following vector function 

g = r'^^u (25) 

is a nontrivial analytic symmetry for the original DS (22) (in the case that p — (x^ + 2/^)^, 
we choose e.g. g — {x ■ y)^u)^ indeed 

{Au + pu, r'^^Iu} = {Au, r'^^Iu} + p{u • Vr'^^)u - r'^^{u • Vp)u = 

and so we can conclude, e.g. from 12) of Theorem 1" and Theorem 1', that the DS can be 
normalized by a convergent transformation. 
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